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Abstract. - We report new experimental results on granular compaction under consecutive 
vertical taps. The evolution of the mean volume fraction and of the mean potential energy of 
a granular packing presents a slow densification until a final steady-state, and is reminiscent to 
usual relaxation in glasses via a stretched exponential law. The intensity of the taps seems to 
rule the characteristic time of the relaxation according to an Arrhenius's type relation . Finally, 
the analysis of the vertical volume fraction profile reveals an almost homogeneous densification 
in the packing. 



Introduction. - Granular matter is a well known example of athermal systems; it means 
systems where classical thermodynamics do not apply since thermal energy (fceT) is insigni- 
ficant compared to the gravitational energy of a macroscopic grain. A static packing of grains 
is therefore in a mctastable state, indefinitely trapped in a local minimum of the total potential 
energy. When submitted to an external perturbation, the system instantaneously acquires an 
extra mechanical energy and then relaxes to a new metastable configuration, which depends 
on the previous one as well as on the nature of the perturbation. This dependence can be 
investigated by implementing at regular intervals identical external excitations on an assembly 
of grains and analyzing the succession of static metastable states explored by the system. This 
is a common experiment in pharmaceutics when compacting powders but it is also a practical 
way to "thermalise" a granular media and to test the conceptual connection between granular 
compaction and the very slow relaxations of out-of-equilibrium thermal systems |f|,|2|. 

The first experiments in this spirit have been carried out in Chicago |5|. Starting 
from a loose packing of beads confined in a tube, a succession of vertical taps of controlled 
acceleration induces a progressive and very slow compaction of the system. This evolution 
is well fitted by the inverse of the logarithm of the number of taps and, after more than 
f0,000 taps, a hypothetic steady state is still not reached. These results have motivated 
many theorical and numerical works, most of them dealing with the notion of free volume and 
geometric constraint ||. Some of them underscore structural aging effects, as currently 
noticed in glasses. 
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In this paper, we present new compaction experiments in what we believe to be more 
general conditions. Indeed, the previous experiences |||],|^] were realized in a thin cylinder of 
diameter D = 1.88 cm filled with monodisperse glass spheres of diameter d = 1, 2, or 3 mm, 
that is to say a horizontal gap of 10 to 20 beads (Nh ~ 10) between the lateral walls. This 
condition allows a local measurement of the volume fraction with a capacitive method and 
prevents any convection in the packing. But, in return, the boundary effects are very strong 
and may be in particular responsible for the highest values of the volume fraction obtained 
in some Chicago experiments Q, significantly above the random close packing limit (approx- 
imately 64%) which corresponds to the maximal volume fraction in a disordered packing of 
identical hard spheres. 




Fig. 1 - (a) Temporal evolution of the mean volume fraction {$} for different tapping intensities F: 
T = 0.96 (A), T = 1.16 (■), T = 1.60 (O), r = 2.75 (•), and T = 5.0 (x). (b) Temporal evolution of 
the mean potential energy (Z) for the same values of V. To make the curves easier to examine, only 
an approximate 25% of the experimental points of the mean volume fraction are plotted. 



Experimental set-up. - Our experimental setup is the following: a glass cylinder of diam- 
eter 10 cm, filled with 1 mm diameter glass beads on about 10 cm height, is shaken at regular 
intervals (At = 1 s) by an electromagnetric exciter delivering vertical taps, each of them con- 
sisting of an entire cycle of a sine wave (frequency / = 30 Hz). The negative peak acceleration 
a-max felt by the whole system is measured by an accelerometer at the bottom of the cylinder 
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so as to parametrize the tap intensity by the dimensionless acceleration T — a max /g. By 
measuring the absorption of a 7-ray beam through the packing, it is possible to estimate the 
average volume fraction in the bulk ($) as well as the vertical density profile providing 
a local analysis of the packing structure. We can also deduce the mean potential energy of 
the heap: (Z) = / °° z$(z)dz / J °° <f>(z)dz. 

To restrict the boundary effects, we use here a large horizontal gap (Nh ~ 100) even 
if, therefore, we allow convection to occur during the compaction of the beads, giving rise 
in particular to an instability of the horizontal free surface Q. In comparison with the 
previous experiments j^, || |E| , another important difference concerns the vertical pressure in 
the static packing: in the narrow tube used in Chicago's set-up, the pressure felt by the 
heap at a given height does not correspond to the total weight of the upper packing; part 
of this weight is screened by the lateral walls ("Janssen effects"). As a consequence, in a 
static configuration, the vertical pressure is homogeneous in almost all the packing whereas, 
here, the vertical pressure is definitely not saturated in the packing but probably close to the 
hydrostatic pressure. Nevertheless, we do not know to what extent the initial static situation 
can play a role in the dynamical process induced by a tap. 

The measure is deduced from the transmission ratio of the horizontal collimated 7 beam 
through the packing: T — A/Aq where A and Aq are respectively the activities counted on 
the detector with and without the presence of the beads in the cylinder. From the Beer- 
Lambert's law for absorption, we can derive an estimation of the volume fraction in the probe 
zone: $ « —(/.iD)~ 1 ln(T). Here /i is the absorption coefficient of the beads; it was evaluated 
experimentally to fj, « 0.188 cmT 1 for our 7 beam of energy 662 keV ( 137 Cs source). 

The collimated 7 beam is nearly cylindrical with a diameter of 10 mm and intercepts 
perpendicularly the vertical axis of the cylinder of beads. An acquisition-time of 60 seconds 
for each measure was found as a good consensus between the intrinsic uncertainty of the 
radioactive beam and the total duration of an experiment. We then achieve a precision 
A<£> w 0.003. Two types of measurement are used: vertical profile $(z) and mean volume 
fraction (<j>). The vertical volume fraction profile is deduced from 63 measures of the 7- 
transmission at successive heights z with a regular step Az w 2 mm; each measure is an 
average on a horizontal slice obtained by 1 turn rotation of the cylinder about its axis during 
the measure. ($) is estimated from the transmission ratio T averaged on approximately 7 
cm-height from the bottom of the cylinder: the cylinder achieves a vertical translation of 7 cm 
combinated to a rotation of 2 turns so as to permit a significant saving of time. With the aim 
of limiting the duration of the experiments and of avoiding redundant information due to the 
very slow evolution of the system, the measurements are spaced out in time (on a logarithmic 
scale) with 2 measures of profile and 50 measures of (<&) per decade (except 10 for the first 
decade). 

Compaction dynamics. - Several compaction experiments were carried out for different 
values of the tapping strength T in the range [0,6]. Part of the results are presented in figure 
which shows the evolution of the mean volume fraction ($) and of the mean potential energy 
(Z) during 10,000 or 100,000 taps and for a few values of T. These curves are either raw 
data or an average on 2 or 3 realizations. Here we call "time", t, the number of taps and the 
"dynamics" is the succession of static equilibrium induced by the taps. (<&) and (Z) are plotted 
versus log(t + 1) as a convenient way to include the initial state (t = 0) on the logarithmic 
axis. This initial state corresponds to a loose packing (($) = 58.3 ± 0.3%) preparated in a 
reproductible way. 

The typical evolution is a slow compaction of the packing charaterised by an increase 
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Fig. 2 - (a) Volume fraction of the steady-state ($) ss versus V. (b) The experimental estimation 
of the relaxation time t c (■) and the characteristic time Tf (O) of the streched exponential fit as 
functions of the inverse of the tapping intensity V (r c and r/ are defined in the text). The lines are 
linear fits corresponding to Arrhenius laws (see equation (0)). 



of ($} and a reciprocal decrease of (Z). Then, after a varying relaxation time, the system 
finally reaches a steady-state. Note that for T = 0.96, the temporal window accessible in the 
experiments becomes too small to observe the whole relaxation process; but a compaction 
obviously occurs although T is smaller than the "dynamical" threshold T* (T* m 1.2) above 
which there is a collective takeoff of the packing from the bottom of the cylinder. When 
reached, the final steady-state is all the more compact (i.e. small value of (Z) ss and high 
value of ($)ss) as the tapping strength is slight (see figure ||a); but, in return, the number of 
taps needed increases significantly. In all our experiences, the volume fraction stays below the 
random close packing limit which tends to reject any hypothesis of ordering or crystallization 
in the packing. Convection seems also to play a role in the compaction process: on figure 
, we observe indeed a significant change in the dependance of ($) ss with T which might 
correspond to different convective regimes. Under a threshold T c w 2, the final state of the free 
surface of the packing is an inclined plane and indicates a spontaneous breaking of symmetry. 
Above T c , the free surface heaps up moderately and finally takes a flat conical shape probably 
brought about by a nearly torical convective roll. These sorts of free surface instabilities have 
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already been observed, see for instance 
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Fig. 3 - Collapse of the compaction curves obtained with V — 1.01, 1.16, 1.33, 1.48, 1.60, 1.76, 1.95, 
2.15, 2.46, 2.75, 3.15, 3.58, 4.20, 5.0, and 6.0: x = ({$)ss - ($)(*))/(($)«, - (*}(0)) is plotted versus 
(t/rj) 13 . The solid line is the function /(it) = exp(~u) corresponding to the stretched exponential 
law. 



In comparison with the previous experimental results ||, some sharp differences appear. 
First, the shapes of the compaction curves differ significantly, especially concerning the long 
time behavior and the obtaining of a final steady-state which is here definitely established and 
may correspond to a dynamical balance between compaction and convection. Moreover, for a 
given intensity T, the dynamics of the compaction seems slower in the Chicago's experiments, 
particularly for the highest values of T. We believe that these differences are principally due to 
the disparity of the lateral constraint (horizontal gap Nh) in the two configurations (Nh ~ 10 
against Nh ~ 100). To extend the comparison, we have tried to fit ($) with the empirical law 
initially proposed in & 



<*>(*) = Woo - 



A($) c 



l + Sln(l + t/r) 



(1) 



The result is satisfactory concerning the beginning of a typical compaction curve but fails 
to correctly fit the final relaxation up to the steady-state. In particular, the parameter ( < I > ) o 
overestimates sharply the steady-state volume fraction ($) ss at small T. On contrary, our 
datas concerning (<&} as well as (Z) are in very good agreement with a stretched exponential 
function (equation^)) on the whole temporal range. First used by Kohlrausch in 1854 1 10|, this 
expression was far later popularized by Williams and Watts pd| ] . It is now frequently applied 
to a large range of relaxations in disordered thermal systems as glasses (see for example |Q 
and references therein) and is often called KWW law. 



X{t) =X OQ - (X^ - X Q ) exp ( - (t/r f f) with X - ($> or (Z) 



(2) 



As Xqo can be approximated by X ss and X by X(t = 0), equation (||) has only two free 
parameters (tj and (3) . Examples of this fit are presented in figure [l] where the solid lines are 
the stretched exponential laws. For F = 0.96 and more generally for r < 1, only the beginning 
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of the relaxation is accessible and there is so no evidence that a stretched exponential law still 
describes the compaction dynamics; consequently the corresponding fit is plotted in dotted 
line. The values obtained for (3 are in the range 0.5 — 0.8 and tend to increase slightly with T 
whereas r/ decreases strongly with T. 

A general collapse of the compaction curves can be obtained by use of the function \i 
defined as the rate of increase of ($}: X (t) = (($} ss - ($)(i))/(($) ss - ($)(0)) = f{it/T f f). 
Figure || presents the plot of \ versus (t/rj)^; the solid line plots the KWW fit. The fluctu- 
ations at large t correspond to the final steady-state and are particularly significant at large 

r. 

Relaxation time. - To quantify more precisely the influence of V on the dynamics of the 
compaction, it is possible to estimate a characteristic relaxation time. To do this, we use 
two different initial packings: the loose one, already presented above (($) = 58.3 ± 0.3%), 
and a more compact one ((<!>) = 63.2 ± 0.2%). Submitted to identical taps, the first packing 
densities whereas the other dilates and both of them progressively meet, reaching the same 
steady-state. From this meeting point, we can evaluate a time of convergence r c . As most 
of the usual disordered packings of monodisperse beads have a packing fraction comprised 
between 58.3% and 63.2%, this time can also be regarded as a memory effect and can be 
interpreted as the number of taps required for a packing to "forget" its initial configuration. 
When analysing the dependance with T of this characteristic time r c as well as of the parameter 
time Tf of the stretched exponential fit, we found that an Arrhenius behavior (equation(|J)) 
describes reasonably well the experimental dynamics as illustrated in figure ||b: 

r Ci/ (r) = r exp(^) (3) 

In the two cases, we have obtained (r « 9.6, t ~ 6.7) for r c and (r « 8.0, r « 0.9) for r/. 




Fig. 4 - Quasi-homogeneous compaction of the vertical volume fraction profile. The profiles are 
represented in the range 56% < <£> < 62% and < z < 9 cm and correspond to the following tap 
numbers: t + 1 = 1, 3, fO, 32, 100, 316, 1.000, 3.162, and 10.000. 
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Finally, we present in figure [| a typical evolution of the vertical volume fraction profile in 
the case of T = 2.15. To observe the progressive densification of the packing, we use a zoom 
on the zone of interest (56% < $ < 62%) and a three-dimensional representation with the 
third axis corresponding to time on a logarithmic scale. The profile is plotted approximately 
every half-decade. 

We note that the compaction is rather homogeneous in the bulk; there is no upward or 
downward densification front. The profile continuously reaches a final steady-state. This 
asymptotic profile is nearly uniform, a slight positive gradient (d^/dz > 0) appears for T > 2; 
this observation agrees with previous experimental || and numerical Q| results. We have also 
verified that the same steady-state profile is obtained when starting from the initial dense 
packing (($) = 63.2 ± 0.2%) instead of the loose one. 

Conclusion and perspectives. - In conclusion, we have shown that, in the case where the 
lateral constraint is weak, granular compaction by vertical taps is quite similar to a typical 
relaxation of a out-of-cquilibrium thermal system. Indeed, the densification curves can be 
reasonably well fitted by a stretched exponential or KWW law. The characteristic time of 
compaction (r c or r/) follows an Arrhenius relation where the dimcnsionless acceleration T 
plays the role of temperature. These results confirm and reinforce the analogy between com- 
paction dynamics and "glassy" phenomena. Moreover, the acquisition of the vertical volume 
fraction profile has established that the compaction is rather homogeneous in all the packing. 
Further investigations are in progress to study more precisely the influence of the lateral con- 
straint (Nh) and to clarify the link between convection and compaction. 



* * * 

We are grateful to A. Valance and R. Delannay for a careful reading of the manuscript, 
and to S. Bourles and P. Chasle for technical assistance. 

REFERENCES 

[I] Edwards S. F. and Oakeshott R. B. S., Physica A, 157 (1989) 1080; Mehta A. and Edwards 
S. F., Physica A, 157 (1989) 1091 

[2] Nicodemi M., Coniglio A. and Herrmann H. J., Phys. Rev. E, 55 (1997) 3962; Kurchan J., 

J. Phys: Cond. Matt, 12 (2000) 6611 
[3] Knight J. B., Fandrich C. G., Lau C. N., Jaeger H. M. and Nagel S. R., Phys. Rev. E, 

51 (1995) 3957 

[4] Nowak E. R., Knight J. B., Povinelli M. L., Jaeger H. M. and Nagel S. R., Powder 
Technology, 94 (1997) 79 

[5] Nowak E. R., Knight J. B., Ben-Naim E., Jaeger H. M. and Nagel S. R., Phys. Rev. E, 
57 (1998) 1971 

[6] Barker G. C. and Mehta A., Phys. Rev. E, 47 (1993) 184; Philippe P. and Bideau D., Phys. 

Rev. E, 63 (2001) 051304 
[7] Barrat A. and Loreto V., J. Phys. A, 33 (2000) 4401 
[8] Krapivsky P. L. and Ben-Naim E., J. Chera. Phys., 100 (1994) 6778 

[9] Evesque P. and Rajchenbach J., Phys. Rev. Lett., 62 (1989) 44; Laroche C, Douady S. 

and Fauve S., J. Phys. (Paris), 50 (1989) 699 
[10] Kohlrausch R., Pogg. Ann. Phys. Chem., 91 (1854) 179 

[II] Williams G. and Watts D. C, Trans. Faraday Soc, 66 (1970) 80 
[12] Phillips J. C, Rep. Prog. Phys., 59 (1996) 1133 



